The pasta phase within density dependent hadronic models 
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In the present paper we investigate the onset of the pasta phase with different parametrisations of 
the density dependent hadronic model and compare the results with one of the usual parametrisation 
of the non-linear Walecka model. The influence of the scalar-isovector virtual S meson is shown. 
At zero temperature two different methods are used, one based on coexistent phases and the other 
on the Thomas-Fermi approximation. At finite temperature only the coexistence phases method is 
used, npe matter with fixed proton fractions and in /3-equilibrium are studied. We compare our 
results with restrictions imposed on the the values of the density and pressure at the inner edge 
of the crust, obtained from observations of the Vela pulsar and recent isospin diffusion data from 
heavy-ion reactions, and with predictions from spinodal calculations. 
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I. INTRODUCTION 

Frustration is a phenomenon characterized by the ex- 
istence of more than one low-energy configuration. The 
pasta phase is a frustrated system [3, Jl, • At densi- 
ties of the order of 0.006 - 0.1 fm"'' [j| in neutral nu- 
clear matter and 0.04 - 0.065 fm~^ Q in /3- equilibrium 
stellar matter, a competition between the strong and 
the electromagnetic interactions takes place leading to 
a frustrated system. The basic shapes of these complex 
structures were first named [2] after well known types 
of cheese and pasta: droplets (bubbles = Swiss cheese) , 
rods = spaghetti (tubes = penne) and slabs (lasagna) for 
three, two and one dimensions respectively. A droplet 
(bubble) and a rod (tube) have densities larger (smaller) 
than their surroundings, and are normally defined within 
a Wigner-Seitz cell. The pasta phase is the ground state 
configuration if its free energy per particle is lower than 
the corresponding to the homogeneous phase at the same 
density. The pasta phase is expected to exist somewhere 
between a solid and a liquid phase, more like a liquid 
crystal Q. Its mechanical and thermal properties are 
likely to depend on its shape and this study still remains 
to be done. 

These pasta shapes at sub-nuclear densities are ex- 
pected to exist both in the crust of neutron stars (zero 
temperature, very low proton fraction, matter in /3- equi- 
librium) and in supernova (finite temperature, proton 
fraction around 0.3). In neutron stars the pasta phase 
coexists with a neutron gas; in supernova there is no 
neutron gas or it is very low in density . 

In a recent work [4] we have studied the existence of 
the pasta phase at zero and finite temperature within 
three different parametrisations of the relativistic non- 
linear Walecka model (NLWM) namely NL3 [9], TMl 
[Tot GM3 the last one generally used in the 

studies of stellar matter. At zero temperature two dif- 



ferent methods were used: the coexisting phases (CP) 
and the Thomas-Fermi (TF) approximation. We have 
checked that while the final equation of state (EoS) ob- 
tained with the different methods do not vary much, the 
internal structure varies considerably. The TF approxi- 
mation was performed to test the much simpler CP cal- 
culation and we have seen that the success of the CP cal- 
culation depends on the parametrisation of the surface 
energy for very small proton fractions and close to the 
transition densities. At finite temperature only the CP 
method was used and compared with predictions from 
spinodal calculations. The pasta phase shrinks with the 
increase of the temperature and we have found that ho- 
mogeneous matter can be the preferential phase also at 
very low densities depending on the temperature and the 
proton fraction. If equilibrium is imposed the pasta 
phase does not appear in a CP calculation. This indi- 
cates the necessity to use a good parametrisation for the 
surface energy which is temperature, proton fraction and 
geometry dependent, as also stressed in [T^ [T3|. 

The authors of [l^l have related the fraction of the mo- 
ment of inertia contained in the crust of the Vela pulsar 
with the mass and the radius of the neutron star and the 
pressure and density at the crust-core interface. From 
realistic EoS they have obtained an expected range of 
values for the pressure at the inner edge of the crust and 
therefore also a relation between the radius and mass of 
the pulsar. This work shows the importance of under- 
standing the exact density limits of the pasta phase and 
its consequences on the choice of appropriate equations 
of state. More recently a new radius-mass relation for 
the Vela pulsar was obtained taking as constraints recent 
isospin diffusion data from heavy- ion reactions Q . In this 
work both the thermodynamical and the more accurate 
and reliable dynamical method were utilized in order to 
constrain the densities and related pressures of the pasta 
phase present in the crust of neutron stars. The pressures 
were obtained from the equation of state for neutron-rich 
nuclear matter constrained by isospin diffusion data ob- 
tained in the same sub-saturation density range as the 
existing ones in the neutron star crust [l5j. 

It is, however not clear how good are the predictions 
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for the transition density obtained from spinodal calcula- 
tions. Clusterisation of the crust may have been formed 
through equilibrium processes and it is important to com- 
pare spinodal results with equilibrium results, obtained 
from the minimization of the free energy. In [4] a first 
comparison was done and it was shown that as a rule 
the transition densities obtained within an equilibrium 
calculation are larger than the ones determined from the 
dynamical spinodals. 

In the present work we use the same approximations 
(CP and TF) used in to obtain the pasta structures, 
but improve on the choice of the relativistic models, 
i.e., we obtain results with various density dependent 
hadronic models and investigate the influence of the delta 
mesons. We next justify our choices. 

Density dependent hadronic models [H, [l3| have 
shown to provide richer and different results in many 
cases as compared with the simpler NLWM parametri- 
sations [H, [1^, [13]. In many situations the results are 
similar to the ones obtained with non-relativistic Skyrme- 
type models [H iH HI H!] . 

The inclusion of the isovector-scalar virtual (5(ap(980)) 
meson in hadronic effective field theories [U, [2^ influ- 
ences the calculation of the effective masses with impor- 
tant consequences on the symmetry energy, spinodals [l^ 
and other quantities possibly related with the appearance 
of the pasta phase. The 5 field introduces in the isovec- 
tor channel the structure already existing in the isoscalar 
channel, i.e., a balance between a scalar (attractive) and 
a vector (repulsive) potential. 

In the following we consider three density dependent 
coupling parametrisations, TW [lB|, DDH(5 [l3| and 
GDFM J^] and two models with constant couplings, NL3 
and NL3(5. Neither NL3 nor TW include the 5 meson. A 
comparison is done between the transition pressures and 
densities from the pasta phase to homogeneous matter 
obtained within the above mentioned models, with the 
predictions obtained in [5^]. 

The paper is organized as follows: in section II we 
briefly review the formalism underlying the homogeneous 
neutral npe matter. In section III the pasta phase is 
built with the help of the coexisting phases method and 
in section IV with the Thomas-Fermi approximation. In 
section V our results are displayed and commented and 
in section VI our conclusions are drawn. 



II. THE FORMALISM 

We consider a system of protons and neutrons with 
mass M interacting with and through an isoscalar-scalar 
field (t> with mass TOs, a isoscalar- vector field with 
mass m^^ , an isovector- vector field b'^ with mass and 
an isovector-scalar field 5 with mass mg. We also include 
a system of electrons with mass rUe- Protons and elec- 
trons interact through the electromagnetic field A^. The 
Lagrangian density reads: 



C=Y^ C^+ Ce+C^+C^+Cp+Cs+C^, (1) 

where the nucleon Lagrangian reads 

A = ^a7M^^^-M*]V',, (2) 

with 

iD^' ^ ia^^ _ r^y^ - • b'^ - e^^A^, (a) 

M* = M -Ts(t)~VsT-5, (4) 
and the electron Lagrangian is given by 

Ce = [7m {^d"" + - ^e] V'e. (5) 

The meson and electromagnetic Lagrangian densities are 

r —-If PA"' 

where Vl^^^ = d^Yu-d^V^ , B^^, = dph^-dyhp^-V p{h^ x 
hu) and — d^Ai, — 9^/1^. The parameters of the 
models are: the nucleon mass M = 939 MeV, four density 
dependent coupling parameters Fg, F^^, Fp and F^- of the 
mesons to the nucleons, the electron mass vrif. and the 
electromagnetic coupling constant e = ^47r/137. In the 
above Lagrangian density t is the isospin operator. 
From de Euler-Lagrange formalism we obtain coupled dif- 
ferential equations for the scalar, vector, isovector-scalar, 
isovector- vector, electromagnetic and nucleon fields. In 
the static case there are no currents and the spatial vec- 
tor components are zero. In a complete description 
of the Thomas-Fermi approximation applied to different 
parametrisations of the NLWM is given. As the differ- 
ences arising from the inclusion of the 5- mesons and 
the use of density dependent couplings are small, we do 
not repeat the equations here. A rearrangement term is 
the landmark of most density dependent hadronic models 
[28l [29j and the simple mean field approximation (MFA) 
is outlined next so that its appearance is better under- 
stood. The equations of motion for the fields can be 
obtained and solved self-consistently in the MFA (the 
photon and meson fields are classical fields), neglecting 
states of negative energy (no-sea approximation) [ig . 

The meson fields within the mean field approximation 
are obtained from the following equations: 

m^0o = ^sPs, (6) 
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(7) with 



r 

2l P 



's3- 



(8) 



(9) 



The second members of the above equations include the 
the equiUbrium densities p = Pp + Pn, P3 = Pp — Pn, Ps = 
Psp + Psn and ps3 — Psp — Psn where the proton/neutron 
densities are given by 

Pi^^ I P^dp{f,+ -/,_), i=p,n (10) 



and the corresponding scalar density by 



M* 



p2 + M*2 



= (/,+ +/,_), (11) 



with the distribution functions given by 



l + exp[(e*(p)T^^«)/r] ' 



(12) 

i\<^i yyj -I- '^n/-^ J 

where e* = -^p^ + Mf, 

M* = M -T, (bo-T3,rs Ss, (13) 
and the effective chemical potentials are 



l^^^P^- r^Vo ~ T3, bo - S^, (14) 

= ±1 is the isospin projection for the protons and 
neutrons respectively. The density dependent models in 
the mean field approximation contain a rearrangement 
term E« [I^: 



_ dT, , dTp bo dVs 



Ps3 S3 



Notice that for T = MeV the distribution function for 
particles given in equation (|12[) becomes the simple step 
function fi = 0{Pp^ — p^) and the distribution function 
for anti-particles vanishes. 

In the description of the equations of state of a system, 
the required quantities are the baryonic density, energy 
density, pressure and free energy. The energy density 
reads: 



(15) 



K^ = ^ J v'dp^p^ + MfiU + U.), (16) 

(17) 
(18) 
(19) 
(20) 

The pressure is given by: 

P= P^ + P.+P^+Ps+Pp, (21) 



£<y = —Vo^ 
£ - 

£ ^b^ 



i—n,p 



with the partial pressures associated with the nucleons 
and the various fields 



P 



p2 + Mf 



ml f p (9r, 



P, = -^6^(1 + 2 



P dV, 



^P dp 



P& = -^Si (1 + 2 



P dVs 

2 ' dp 

The free energy density is defined as: 

T = £-TS, 
with the entropy density : 

S = ^{£ + P - PpPp - flnPn)- 



(22) 



(23) 



As for the electrons, their density and distribution 
functions read: 



with 



/e±(r,p,t) = 



1 + exp[(ee T Ate)/T] ' 



(24) 



(25) 



i—n,p 



where /le is the electron chemical potential and ee = 
■\/p^ + mf . We always consider neutral matter and there- 
fore the electron density is equal to the proton density. 
In the calculation of the non-homogeneous phase the 
Coulomb energy of the proton and electron distributions 
is included. We study both matters (homogeneous and 
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pasta structured) , with a fixed proton fraction, as we get 
in fieavy ion collisions (although in this case matter is 
not neutral) and in /3-equilibrium as in stellar matter. In 
the latter, charge neutrality conditions fix the electron 
chemical potential and density. The onset of muons oc- 
curs above the transition density to homogeneous phase 
and therefore the proton density remains equal to the 
electron density. 

The energy density, pressure, free energy density and 
entropy density of the electrons are 



+ mUfe+ + fe-), (26) 



(/e++/e-), (27) 



(28) 



bulk properties, using the following parametrisation for 
the a and lu mesons: 



with 



and 



Ti{p) = Ti{po)hi{x), x^p/po 



lli[X) ^ ttiZ ; TTTT, l = S,V, 



^l + c,(x + d,)'' 



hp{x) — exp[— ap(x — 1) 



(30) 



(31) 



(32) 



and 



for the TW model In the case of the DDHJ model 
[l3| we use the TW parametrisation for To- and and 
for the p and 6 mesons we take: 

hi{x) = aiexp[-6i(x - 1)] - Ci{x - di), i ^ p, 6. 

The parameters ai,bi, Ci and di are given in Table El and 
Po is the saturation density. 



Pe - PePe)- 



(29) 



To obtain the equations for the TW parametrisation 
[l^ of the density dependent hadronic model, all infor- 
mation on the 5 meson is excluded. For the NL3 [9] 
parametrisation the density dependent parameters are 
substituted by the usual coupling constants £<t, g^j, gp 
and non- linear parameters are included (see for in- 
stance). The NL3(5 parametrisation is defined with the 
same values for g„- and g^j as in the NL3 parametrisa- 
tion, gp = 14.29 and gg = 7.85, in such a way that the 
symmetry energy has the same value at the saturation 
density as the NL3 parametrisation. We show in Table |T] 
the nuclear matter properties reproduced by the models 
we discuss in the present work. 



TABLE I; Nuclear matter properties at the saturation den- 
sity. 





NL3/NL35 


TW 


DDH5 GDFM 




[9] 


[16] 


[17] 


[12] 


B/A (MeV) 


16.3 


16.3 


16.3 


16.25 


po (fm-^) 


0.148 


0.153 


0.153 


0.178 


K (MeV) 


272 


240 


240 


337 


S.ym. (MeV) 


37.4 


32.0 


25.1 


32.11 


M'/M 


0.60 


0.56 


0.56 


0.68 


L (MeV) 


117/148 


55 


44 


57 


n (MeV) 


14.55 


15.18 


15.18 


15.95 



The density-dependent coupling parameters are ad- 
justed in order to reproduce some of the nuclear matter 



TABLE II: Parameters of the TW \W\ and DDH5 models 

02.. OS] 



i 


m,(MeV) 


r»(po) 


ai 


bi Ci 


di 


a 


550 


10.72854 


1.365469 


0.226061 0.409704 


0.901995 


uo 


783 


13.29015 


1.402488 


0.172577 0.344293 


0.983955 


pTW 


763 


7.32196 


0.515 






Pddhs 


763 


11.727 


0.095268 


2.171 0.05336 


17.8431 


5 


980 


7.58963 


0.01984 


3.4732 -0.0908 


-9.811 



The density dependent parametrisation GDFM ob- 
tained in [l^ takes into account the renormalization of 
the relativistic mean-field theory due to Fock exchange 
terms. It ensures a good description of the properties of 
the equation of state at high density as obtained with cal- 
culations for asymmetric nuclear matter [isj with Dirac- 
Brueckner-Hartree-Fock calculations. 

The GDFM parametrisation for all four mesons cou- 
pling parameters reads: 



ri{p) = ai + {bi + dtx^) exp(-Cix). 



(33) 



Around the saturation density a correction to the cou- 
pling parameter for the lo meson is introduced: 



t{p) = '['ujip) - flcor exp 



P- PO 



(34) 



where Ocor — 0.014 and 6cor = 0.035 fm^'^. 

The properties we discuss in the following depend on 
the isovector channel of the nuclear force, mainly the re- 
sults obtained for /3-equilibrium matter. Therefore we 
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TABLE III: Parameters of the GDFM model hj 



III. COEXISTING PHASES 



i 


m,(MeV) 






Ci 


di 


a 


550 


7.7868 


2.58637 


2.32431 


3.11504 


Ul 


782.6 


9.73684 


2.26377 


7.05897 




P 


769 


4.56919 


5.45085 


1.20926 




S 


983 


2.68849 


6.7193 


0.503759 


0.403927 
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FIG. 1; a) Symmetry energy and b) slope of symmetry energy 
for the models under study. 



show in Fig. [T] the symmetry energy and the slope of 
the symmetry energy L = 3 pa desym/dp. This quantity 
defined at the saturation density is given in Table |T1 The 
symmetry energy of NLW models becomes quite hard for 
densities above ^ 0.1 fm^"^. However, at subsaturation 
densities, the 5 meson gives rise to a softer symmetry 
energies: this is true both for NL3(5 and for DDH(5 and 
GDFM, DDH(5 being softer than GDFM. Looking at the 
slope of the symmetry energy we see that GDFM and 
DDH(5 show a very similar behavior. Except for TW all 
models show at low densities a decrease of the slope fol- 
lowed by an increase of the slope above p ^ 0.05 fm^"^ 
for NLW models or p 0.12 fm~'^ for density depen- 
dent models with the 5 meson. For TW the slope always 
decreases with density more slowly than all the other 
models. 



A. Nuclear pasta 

As in I, g, for a given total density p and proton 
fraction Yp = Pp/ p the pasta structures are built with 
different geometrical forms in a background nucleon gas. 
This is achieved by calculating from the Gibbs' conditions 
the density and the proton fraction of the pasta and of 
the background gas, so that in the whole we have to solve 
simultaneously the following eight equations: 



2 ll I 



ml4' r,pi 



mjSi'^Tspii, 



II 



fp'p + (1 - f )Pp =Pp= Ypp, 



(35) 
(36) 

(37) 

(38) 

(39) 

(40) 

(41) 



where I and II label each of the phases, / is the volume 
fraction of phase I: 



P 



II 



nil 



(42) 



and Yp is the global proton fraction. The density of elec- 
trons is uniform and taken as pe = Ypp. For the NL3 
and NL3(5 parametrizations, non- linear a terms must be 
included in ([37]) and ([38 |l . 

The total pressure is given hy P = + Pe- The total 
energy density of the system is given by 



£ = f£' + f)£'' +£e+ Esurf + 



Coul 1 



(43) 



where, by minimizing the sum Egurf +£coui with respect 
to the size of the droplet/bubble, rod/tube or slab we get 
3] Ssurf = "^Scoui, and 



2a 



//^^2/3 



Seoul = ^(e27r<i>)i/3 {aD{p'p - p^^)) 



(44) 



where a — f for droplets and a = 1 — f for bubbles, a is 
the surface energy coefficient, D is the dimension of the 
system. For droplets, rods and slabs. 



f-l-ln(f) 
D+2 ' 



1 

D+2 ' 



^ = 1,3; 



D = 2. 



(45) 
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and for bubbles and tubes the above expressions are valid 
with / replaced by 1 — /. 

Concerning the surface energy, the authors of Q state 
that, in this case, the appearance of the pasta phase es- 
sentially depends on the value of the surface tension. We 
have fixed the surface tension at different values and con- 
firmed their claim. We have parameterized the surface 
energy coefficient in terms of the proton fraction accord- 
ing to the functional proposed in [3ll | , which was obtained 
by fitting Thomas- Fermi and Hartree-Fock numerical val- 
ues with a Skyrme force. 



with 



16 + 6 
(T3i^ + f 



T 



(46) 



(47) 



CTo = 1.03 MeV/fm2 and 6 = 24.4 and is the critical 
temperature above which there is a smooth transition 
from the gas phase to the liquid phase 19] and given 
in Table [J We have checked that small variations of 
this temperature do not affect our results. The proton 
fraction considered throughout the calculation of a is the 
one of the denser phase. 

Each structure is considered to be in the center of 
a charge neutral Wigner-Seitz cell constituted by neu- 
trons, protons and leptons [13] ■ The Wigner-Seitz cell is 
a sphere/cilinder/slab whose volume is the same as the 
unit BCC cell. In [s^ the internal structures are associ- 
ated with heavy nuclei. Hence, the radius of the droplet 
(rod, slab) and of the Wigner-Seitz cell are respectively 
given by: 



R 



D 



ctD 



47re2(p^ 



1/3 



Rw 



Rl 



(48) 



together with ([35]), dST]), ^ and (gOl). Here the 

density of electrons is no longer taken uniform as in the 
last section, but appears as the solution of equation ((52|l . 
The densities of interest to the study of the pasta phase 
are too low for the muons to appear, which generally 
occurs for densities above 0.1 fm""^ |23l ]. 



IV. PASTA-PHASE WITHIN THE 
THOMAS-FERMI APPROXIMATION 

In the present work we repeat the same numerical pre- 
scription given in Q where, within the Thomas-Fermi ap- 
proximation of the non-uniform npe matter, the fields are 
assumed to vary slowly so that the baryons can be treated 
as moving in locally constant fields at each point. In the 
Thomas-Fcrmi approximation, the energy is a functional 
of the density given by: 



Etf = J d^r j ^»(^) 



\[{m,{v)f + mlbl{v)]-\[{V5,{v)f + mj5l{v)] (53) 



r^^o(r)p+-r^6o(r)p3 



(54) 



-[VAo(r)]2 + e(pp-pe)v4o(r) , (55) 



B. Stellar pasta 



where 

I fPFi (r) 

Er^— dpp\p^ + M*Y^^ ,t=p,n , (56) 

Jo 



In this case, hadronic matter is in /3 equilibrium. The 
condition of /? equilibrium in a system of protons, neu- 
trons, electrons and muons is 



— pn Me 7 



(49) 



where = Mm- muons are added, the imposition 

of charge neutrality requires that 



Pp = Pe+ Pti- 

The Gibbs conditions to be enforced are 



Mri Mn 7 Me Me i 



(50) 
(51) 



and 



and 



1 /-PF^Cr) 

= 4 / dppHp' + ml)'/' 

Jo 



(57) 



The definition of the thermodynamic potential reduces 



to 



n = ETF[p^]- J2 ^^^Jdrp^{r). (58) 



i—n,p,e 



The minimization of the above functional with the con- 
straint of a fixed number of protons, neutrons and elec- 
trons yields the equations: 



/ {pI ~Pi- Pi) + (1 - /) {P'Z pV Pi') = (52) bl(r) + Mfir)^/' + r.Vb(r) + '-TMr) + 
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+ e^o(r) = fJ.p 



(59) 



where Af* and M* are given in eq. p^ and 
(p2 (r)+m2)i/2_e^,(r)=^,. 



(60) 



(61) 



The numerical algorithm for the description of the neu- 
tral npe matter was discussed in detail in '4!| . The Poisson 
equation is always solved by using the appropriate Green 
function according to the spatial dimension of interest 
and the Klein-Gordon equations are solved by expanding 
the meson fields in a harmonic oscillator basis with one, 
two or three dimensions based on the method proposed 
in 133. 



V. RESULTS AND DISCUSSIONS 

We show and discuss the results obtained using the co- 
existence phases (CP) and the Thomas-Fermi (TF) meth- 
ods in the framework of the several relativistic models 
presented, always for npe matter. We start with the re- 
sults at T = MeV. 
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FIG. 2: npe matter (pasta -I- homogeneous phase) energy 
per particle at T = MeV for RMF models with constant 
couplings [a), c)] and density dependent coupling models [b), 
d)]: Vv = 0.5 (top), yp = 0.3 (bottom). 

In Fig. [2]results for a homogeneous description of mat- 
ter are compared with the output of the CP calculation 
for two proton fractions. We see that all the models 
predict the existence of a non-homogeneous pasta phase: 
in models with constant couplings, NL3 and NL3(5, this 
phase clearly decreases if j/p decreases (higher asymme- 
tries); a different behavior occurs for density dependent 



coupling models where the non-homogeneous phase ex- 
tends to higher densities and may even increase when 
the proton fraction reduces. From Table IIVI one can see 
that the TF method also predicts for density dependent 
coupling models either a small increase (DDH(5) or just a 
small decrease (TW) of the extension of the pasta phase, 
if we decrease the proton fraction from 0.5 to 0.3. The 
GDFM model presents a large pasta phase at both asym- 
metries but it is with the DDH(5 model that most variety 
of pasta structures can be observed (see Fig. E]). 

We also point out that the models NL3 and NL3i5 as 
well as TW and DDHi5 show the same transition densities 
for z/p = 0.5 because the effect of the S meson is only seen 
for asymmetric matter for the CP and spinodal calcula- 
tions. Both couples of models have the same coupling 
constants for the a and the tu mesons and only differ in 
the isovector channel, namely the p and S meson cou- 
plings. However, in the TF calculation the distribution 
of protons and neutrons is free to adjust itself to the low- 
est energy configuration. As a consequence the proton 
and neutron density distributions do not coincide within 
the Wigner Seitz cell and the p and S fields are not zero. 
However, the differences are not large enough to change 
the transition density and we still get the same transition 
density within TF, for Up = 0.5 and the couples (NL3, 
NL3,5) and (TW, DDH(5). 

For matter in beta-equilibrium at T = MeV the en- 
ergy per particle for the pasta is always slightly larger 
than the corresponding homogeneous matter in the mod- 
els with constant couplings within the CP method; so, in 
these cases, the pasta is never preferred to represent the 
ground state of the system. The absence of a pasta phase 
in the NL3 and NL3(5 parametrisations is related to the 
very high values of the surface tension coefficient a for 
these models as can be seen in Fig. 

•- For the DDH models we show in Fig. [5] that these 
models predict non-homogeneous phases occurring for a 
relatively small range of densities (below p ~ 0.027 fm~'^) 
for the GDFM model but extending up to p 0.09 fm"'^ 
in DDH(5. TW shows an intermediate behavior. Due to 
the low symmetry energy of the DDHJ model, which is 
only 25 MeV at saturation density, it was shown in [2^ 
that /3-equilibrium matter would even present at low den- 

' sities a range of densities with a negative compressibility. 
This behavior is favoring the appearance of a larger va- 
riety of cluster forms. 

Surface tension depends on the proton fraction of the 
high density region inside the Wigner-Seitz cell. Models 
with a larger symmetry energy give a dense region with 
a larger proton fraction and therefore a larger surface en- 
ergy (see eq. p6)) ). We also confirm that DDH(5 shows 
the smallest values for this quantity. Within the TF 
method, on the other hand, the pasta phase is also found 
with the models bearing constant couplings. This is due 
to the fact that in the TF approach the surface energy 
is calculated self-consistently and is not introduced by 
hand. We know, however, that the TF approach predicts 
a too steep surface and, therefore, we may expect that 
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TABLE IV: Transition densities in fm ^ and corresponding pressures (CP and TF calculations) for the non-homogeneous to 
homogeneous phase at the inner edge of the crust (T = 0)MeV. 



model 


EoS 


dynamical thermodynamical 


pasta (CP) 


P (CPl 


pasta (TF) 


P fTF") 






bJJllilJLlcll 


c 1 "Ki /H a \ 
bJJlliULlcll 




MeV/fm^ 




MeV /fm^ 










uniform matter 




uniform matter 




NL3 


V-r. — 5 


U.Uoo 


u.uyo 


0.096 


2.65 


0.112 


2.64 


NL3 


7;„ — n 3 


U.UoU 


u.uy'4: 


0.079 


1.05 


0.100 


1.14 


NL3 


/?-equil. 




u.uuo 






0.054 


0.24 


NL35 


V-r. — 5 

yp 


U.Uoo 


u.uyo 


0.096 


2.65 


0.112 


2.64 


NL3(5 


?;„ — 3 

yp 


n n7Q 

u.u i y 


u.uyo 


0.080 


1.06 


0.099 


1.14 


NL3(5 


/S-equil. 


OAS 


U.UOU 






0.051 


0.16 


TW 


yp ^■'J 


U.UOO 


u.uyu 


0.098 


2.74 


0.113 


2.66 


TW 


tVn = 3 


0.084 


0.095 


0.099 


1.40 


0.111 


1.27 


TW 


/3-equil. 


0.075 


0.085 


0.060 


0.26 


0.076 


0.40 


DDH(5 


yp = 0.5 


0.083 


0.096 


0.098 


2.74 


0.113 


2.66 


DDm 


yp = 0.3 


0.084 


0.094 


0.107 


1.55 


0.115 


1.29 


DDm 


/3-equil. 


0.079 


0.085 


0.089 


0.17 


0.079 


0.10 


GDFM 


yp = 0.5 


0.133 


0.141 


0.119 


3.55 


0.144 


3.81 


GDFM 


yp = 0.3 


0.131 


0.138 


0.119 


1.79 


0.140 


1.83 


GDFM 


/9-equil. 


0.051 


0.058 


0.027 


0.04 


0.052 


0.13 
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FIG. 3: Comparison of the phase diagrams at T = MeV for a) yp — 0.5 and b) yp = 0.3 obtained with the coexisting phases 
method (CP) for several models. From bottom to top the colors represent droplets, rods, slabs, tubes and homogeneous matter. 



a ^ uantal approach would predict a larger pasta phase 

The spinodal surface gives information about the min- 
imal dimension of the pasta phase. A spinodal decompo- 
sition would be expected in the case of a fast transition, 
however in stellar matter we may expect that there is 
always plenty of time to get equilibration. In [Bj it was 
shown that the thermodynamical spinodal results for pn 
matter did not differ very much from the dynamic spin- 
odal ones for npe matter. This seems to indicate that the 
Coulomb interaction and surface tension do not influence 
a lot the pasta phase extension. The thermodynamical 
spinodal for npe matter either does not exist for den- 
sity dependent hadronic models or is very small for NLW 



models due to the large incompressibility of electrons. 
However, although thermodynamically stable, npe mat- 
ter clusterizes as soon as it suffers a density fluctuation 
due to any kind of perturbation. Therefore we discuss 
the extension of the non-homogeneous phase by analyz- 
ing the dynamical spinodal for npe matter, within the 
several models considered. 

In Fig. [5] we display both the thermodynamical spin- 
odals for np matter and the dynamical spinodals for npe 
matter for a momentum transfer fc = 80 MeV, which es- 
sentially defines the envelope of the spinodal surfaces for 
all k values [13] ■ The dynamical spinodals are smaller 
than the thermodynamical ones as expected. The cross- 
ing density of the EoS with a fixed proton fraction equal 
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FIG. 4: Surface tension coefficient cr at T = MeV calculated 
according to eq. (I46p . 
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FIG. 5: npe matter energy per particle at T = MeV for DDH 
models at /3 equilibrium. The full, dashed and dot-dot dashed 
lines stand for homogeneous matter and the corresponding 
dotted lines for pasta phases. 
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FIG. 6: a) Thermodynamical spinodal for np matter and b) 
dynamical spinodal for npe matter, for a transfer momentum 
fc = 80 MeV, for all the models under study. 



to 0.5 and 0.3 or for /3-equilibrium matter is given in 
Table HVl 



A larger extension of the pasta phase within the 
GDFM model is expected from the thermodynamical 
spinodal which we show in Fig. [S) We notice that GDFM 
has a very peculiar behavior with a much larger thermo- 
dynamical spinodal. There is, however, an intermediate 
density region where matter is not so unstable and the 
presence of electrons is enough to raise the instability giv- 
ing origin to two disconnected unstable regions. Compar- 
ing all the spinodals we expect smaller non-homogeneous 
regions for NL3(5 and a larger one for GDFM if yp is not 
too small. For very asymmetric matter like matter in /3- 
equilibrium, the TW and DDH(5 models bear the largest 
pasta phases and NL3, NL3(5 and GDFM predict simi- 
lar results. We also verify that the dynamical spinodal 
predicts a slight increase of the unstable region when j/p 
decreases from 0.5 to 0.3. This behavior is directly re- 
lated with the concavity of the spinodal at yp = 0.5. In 
[23 | it was seen that the concavity of the thermodynam- 
ical spinodal for the TW parametrization at yp = 0.5 
is smaller than the one obtained with the NLWM. The 



presence of electrons and the Coulomb field in the calcu- 
lation of the dynamical spinodal gives rise to a spinodal 
that is not symmetric with respect to the yp = 0.5 axis. 
The spinodal may extend to larger densities for smaller 
proton fractions and the same isospin asymmetry. As 
discussed in we also expect a larger extension of the 
non-homogeneous phase if the electron contribution is 
described correctly, and the Coulomb field included self 
consistently, what stabilizes npe matter and extends the 
non-homogeneous phase. 

It is interesting to compare the density transitions ob- 
tained within the spinodal approaches with the corre- 
sponding values determined from the minimization of the 
free energy both within the CP and TF approaches. As 
discussed before, the dynamical spinodals are expected 
to indicate a lower limit. Within the present models, the 
transition densities obtained from the thermodynamical 
spinodal are ~ 10-15% larger than the values obtained 
from the dynamical spinodal, similar to the conclusion 
drawn in Q. For the proton fractions 0.5 and 0.3 these 
values are always larger than the ones obtained within 
an equilibrium calculation, either CP or TF, except for 
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FIG. 7: Phase diagrams at T = MeV and t/p = 0.5 obtained 
with the TF method for several models. From bottom to top 
the colors represent droplets, rods, slabs, tubes, bubbles and 
homogeneous matter. 



the GDFM model for which the CP results are smaller 
than the spinodal ones. For the /3-equilibrium calculation 
the CP method predicts no pasta phase for the NLWM 
parametrisations (NL3 and NL35) . This is due to the non 
self-consistent description of the surface in the CP ap- 
proach. The TF approach, which treats self-consistently 
the surface, predicts, for all models, a transition den- 
sity larger, but very similar, to the one predicted by the 
dynamical spinodal calculation. This result is very inter- 
esting because it implies that to calculate the transition 
density at the inner edge of the compact star crust it is 
enough to use a dynamical spinodal calculation. 

In Table|V]we compare the transition densities between 
the different pasta geometries obtained in the present cal- 
culation with the results from jl^] , where both a TF and 
a microscopic calculation were done. In the TF calcula- 
tion the surface description was not fully self-consistent 
because it involved the inclusion of a surface energy pa- 
rameter that was adjusted to reproduce the experimen- 
tal binding energy of the nucleus ^°^Pb. We conclude 
that with a self-consistent TF calculation the transition 
densities between the different geometries and from non- 
homogeneous to homogeneous matter are quite smaller 
than the results in [l3| obtained within a microscopic de- 
scription of the pasta structures including pairing effects. 
This comparison should also be done for different proton 
fractions and not only for /3-equilibrium matter. 

^From the phase diagrams shown in Figs. [3^) and[3}3) 
we observe that most models predict the formation of 
inhomogeneities of the type droplet, rod, slab and tube 



TABLE V: Transition densities in fm between the different 
geometries at T = MeV and for the GDFM model. Gom- 
parison with the results taken from 12]. 





GP 


TF 


TF [12] H [12] 


droplet-rod 




0.047 


0.048 0.052 


rod-slab 




0.048 




slab-hom. 


0.027 


0.052 


0.061 0.064 



for the asymmetries considered. For Up = 0.3 only NL3 
evolves to homogeneous matter without a tube- like struc- 
ture and, on the other side, DDH(5 predicts the appear- 
ance of bubbles in a narrow interval of densities. These 
differences are due to the dependence of the surface en- 
ergy on the proton asymmetry and on the slope of the 
symmetry energy. In Fig. [7] the phase diagrams obtained 
with the TF method for the models under study are dis- 
played with Up = 0.5. In this case the bubble structure, 
not present in Fig. [3^), appears. The transition densi- 
ties are systematically higher with TF than with CP, as 
seen in Table ITVl and, therefore another phase structure 
is accommodated. 

At this point we compare our results shown in Table 
IIVI for npe matter in /3-equilibrium with the predictions 
given in According to ^ the transition density and 
pressure from the liquid core to the solid crust at the 
inner edge of neutron stars should lie within the ranges 

0.04 < p< 0.065 fm"^ 0.01 <P< 0.26 MeV/fm^ 

(62) 

The values for the pressure were constrained to the val- 
ues for the slope of the symmetry energy at the sat- 
uration density given by L = 86 ± 25 MeV. The lim- 
its on the transition pressure defined in [5| a re, how- 
ever, quite smaller than the ones given in [14| . namely 
0.25 <P< 0.65MeV/fm^ obtained from realistic EoS. 

One can see in Table H] that the NL3 and NL3(5 models 
have L values respectively slightly larger and larger than 
the upper limit of the proposed L values. However, using 
the TF method, we get, for the transition densities and 
related pressures from the pasta phase to the homoge- 
neous phase, values within the proposed range. On the 
other hand, the density dependent models have L values 
slightly smaller than the lower limit of the above L range. 

The results for the transition pressure and density in 
the TW model lie just at the border of the rectangle that 
sets the limits given in eqs. for the CP calculation 
but they become too high when the TF method is used. 
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FIG. 8: Density and pressure of the inner edge that divides 
the liquid core from the solid crust of neutron stars according 
to d. 
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For the DDH(5 parameter set, the results for the pressure 
are inside the range shown in eas. lp^ both in the CP and 
TF approaches. However, for this model the transition 
density is too high for both calculations. For the GDFM 
parametrisation the density within the CP method is not 
good since it comes up too low, but the TF results lie 
within the imposed constraints. For the NL3 and NL3(5 
parametrisations, the pasta phase is only obtained within 
the TF results and they come out inside the constrained 
range. These observations are summarized in Fig. [51 

Within the dynamical spinodal calculation, the tran- 
sition densities for TW and DDH(5 are high, but all the 
other models are within the density range given in eqs. 

Also from Table [TVl one observes that the influence of 
the 6 meson is only effective for large proton asymme- 
tries. For proton fractions 0.5 and 0.3, NL3 and NL3(5 
give similar results. The effect of the 5 meson is only ob- 
served in the /3-equilibrium matter results: for the con- 
stant coupling models the inclusion of the 5-meson makes 
the pasta phase range a bit smaller, the results of the TF 
calculation being in good agreement with the dynami- 
cal spinodal ones. Among the density dependent models 
and considering all type of calculations presented, we see 
that the extension of the pasta phase for DDH(5 is larger 
than the corresponding one within TW which does not 
include the 6 mesons. This is due to the low value of the 
symmetry energy within the DDH(5 model. 




FIG. 9: npe matter (pasta + homogeneous phases) free energy 
per particle at T = 5 MeV for RMF models with constant 
couplings: NL3 - a) and c), NL3J - b) and d). j/p = 0.5 (top). 
Up — 0.3 (bottom). 

We now comment on the results with finite tempera- 
ture, all of them obtained with the CP method. In Fig. 
[9] we plot the free energy per particle for the models with 
constant couplings. As expected from the T — MeV 
results, no pasta phase appears when /3-equilibrium is 
enforced. For fixed proton fractions (0.5 and 0.3), the 



pasta phase shrinks with temperature. At very low den- 
sities the homogeneous phase has a lower free energy than 
the pasta phase. 

This behavior had already been noticed in ^4] and is re- 
produced with DDH parametrisations, as seen in Fig. [10] 
for Up — 0.5, 0.3 and for matter in /3-equilibrium. From 
Fig. [To] and [Tl] it is seen that the size of the pasta phase 
depends on the asymmetry of the npe matter and on the 
chosen parametrisation. Just two models provide pasta 
phase within CP for matter in /3-equilibrium at T = 5 
MeV: TW and DDHi5, the second being larger than the 
first. 

In Table IVII the transition densities obtained within 
the CP approach and using the thermodynamical spin- 
odal are given. No data for a dynamical spinodal calcula- 
tion at finite temperature is available, except for NL3, see 
4], where the transition densities 0.080, 0.077 fm^'^ were 
given respectively for yp = 0.5 and 0.3 at T = 5 MeV. 
In this case no crossing occurs for /3-equilibrium matter. 
These values are similar to the ones obtained within the 
CP calculation. As discussed before we expect that the 
dynamical spinodal defines a lower limit for the transition 
density. From the discussion of the results obtained for 
the pasta phase at T = MeV, we also expect that the 
CP calculation only gives a lower limit for the transition 
density because of the description of the surface which 
takes too large values. The thermodynamical spinodal 
only suggests an order of magnitude, which, according 
to [j] and for NL3, was always a bit larger than the val- 
ues coming from the CP approach and closer to the TF 



TABLE VI: Transition densities in fm~"^ and corresponding 
pressures (CP calculation) for the non-homogeneous to homo- 
geneous phase at the inner edge of the crust (T = 5 MeV). 



model EoS thermodynamical pasta (CP) P (CP) 

spinodal versus MeV/fm^ 

versus EoS EoS uniform matter 

NL3 yp = 0.5 0.094 0.089 2.42 

NL3 yp = 0.3 0.090 0.075 1.02 

NL3 /3-equiI. 

NL3(5 yp = 0.5 0.094 0.089 2.41 

NL35 yp = 0.3 0.090 0.072 0.96 

NL3S /3-equil. 

TW yp = 0.5 0.095 0.094 2.59 

TW yp = 0.3 0.094 0.094 0.94 

TW /3-equiI. 0.051 0.035 0.15 

DDm yp = 0.5 0.095 0.094 2.59 

DDH<5 yp = 0.3 0.093 0.103 1.51 

DDH(5 /3-equiI. 0.073 0.042 0.10 

GDFM yp = 0.5 0.140 0.101 2.84 

GDFM yp = 0.3 0.137 0.086 1.20 
GDFM /3-equiI. 0.029 
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FIG. 10: npe matter (pasta + homogeneous phase) free energy per particle at T = 5 MeV for RMF models with density 
dependent couplings: TW - a), d) and g), DDH5 - b), e) and h), GDFM - c) and f). Up — 0.5 (top), j/p = 0.3 (middle) and 
/3-equilibrium (bottom) . 



TABLE VII: Highest temperatures for which pasta phase was 
found in MeV. 



NL3 NL35 TW DDH5 GDFM 

[9] [16] [12] 

j/p = 0.5 12.0 12.0 14.0 14.0 10.0 

yp = 0.3 10.0 7.0 13.0 12.0 8.0 

/3-equil. - - 5.0 8.0 4.0 



results for the NLW models. 

For the DDH models, it is seen that within the CP 
approach, and similarly to the result already discussed 
at T = MeV, the pasta phase at T = 5 MeV does 
not decrease when going from j/p — 0.5 to t/p = 0.3 for 



both the TW and DDH(5 models. It is also seen that 
for symmetric matter the prediction obtained from the 
thermodynamical spinodal is generally quite larger than 
the values obtained within the CP calculation, except 
for the DDH(5 model, when they are similar. This may 
indicate that within a TF calculation larger transition 
densities would be obtained. 

Self-consistent calculations at finite temperature, both 
for the pasta phase and for the dynamical spinodal still 
have to be implemented. These are of particular interest 
because neutrino trapping occurs at finite temperature, 
and we expect that the interaction of neutrinos with the 
npe clusters may affect the neutrino energy deposition in 
stellar matter. However, while the CP calculation fails to 
predict the transition densities for /3-equilibrium matter, 
we expect that it gives reasonable results for the pro- 
ton fractions of interest for stellar matter with trapped 



13 



(a) 



(b) 



E 

Q. 



0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 




NL3 NL3S TW 



DDH5 GDFM 



Horn. 
Matter 



Slabs 



— Rods 

— Droplets 

Horn, 
f Matter 





0.14 




0.12 




0.1 


CO 






0.08 






Q. 


0.06 




0.04 




0.02 




- 



NL3 NL38 TW DDH8 GDFM 



Tu 
bes 



Horn. 
Matter 



Slabs 



Rods 

Droplets 

Horn, 
f Matter 



FIG. 11: Comparison of the phase diagrams at T = 5 MeV for a) yp — 0.5 and b) yp = 0.3 obtained with the coexisting phases 
method (CP) for several models. From bottom to top the colors represent homogeneous (for some cases only), droplets, rods, 
slabs, tubes and homogeneous matter. 



neutrinos, yp ^ 0.3. 

We show in Table IVIII the highest temperatures for 
the existence of the pasta phase for each of the mod- 
els discussed in this work within the CP approach. We 
take these values as lower bounds for the limit tempera- 
ture above which the non-homogeneous phase disappears. 
It is also worth mentioning that we have searched for 
the pasta structures at temperature steps of 1 MeV, due 
to the uncertainties mentioned above, and therefore the 
limiting temperature given is not more than an order of 
magnitude. There are two different situations that had 
to be considered when the homogeneous phase was taken 
as the equilibrium configuration: a) the pasta phase ex- 
ists but it is not the preferential state of matter because 
its free energy comes out higher than the homogeneous 
phase and b) the pasta phase does not exist within the 
precision of our calculations. 



VI. CONCLUSIONS 

In the present work we have investigated the exten- 
sion of the pasta phase for npe matter described within 
rclativistic density dependent models, namely TW [l6[, 
DDH5 [13 and GDFM both at zero and finite tem- 
peratures . The pasta phase was calculated at zero tem- 
perature within Thomas-Fermi and compared with re- 
sults obtained in a simplified approach, the coexistence 
phases (CP) Due to the approximate way the sur- 
face is described within the last approach the pasta phase 
comes out smaller with CP than with TF: as in [j] we 
conclude that a correct description of the surface energy 
and its dependence on the isospin, temperature and ge- 
ometry, is essential to obtain better results using the CP 
formalism. 

The effect of including the (5-meson was also explicitly 
investigated: together with DDH(5 and GDFM we have 
also considered NL3(5. It was seen that models with the 
same description of the isoscalar channel and the same 



symmetry energy at the saturation density, namely the 
couple (NL3, NL35), showed a smaller non- homogeneous 
phase for asymmetric matter when the S meson was in- 
cluded. This effect does not occur for the couple (TW, 
DDH(5) because although both have the same descrip- 
tion of the isoscalar channel, the symmetry energy of the 
DDH(5 at saturation is smaller. As a result, the extension 
of the non-homogeneous phase within the DDHi5 model 
is the largest one for /3-cquilibrium stellar matter and 
is larger than the corresponding non-homogeneous phase 
within the TW model. 

Results were compared with previous studies done 
within NLWM and the predictions obtained from the 
analysis of the thermodynamical and dynamical spin- 
odals. One of the main conclusions is that density de- 
pendent hadronic models generally predict larger non- 
homogeneous phases for asymmetric matter than NLW 
models. In fact, for /3-equilibrium matter a similar con- 
clusion had been taken in [23] only from the analysis of 
the crossing of the /3-equilibrium EoS with the dynamical 
spinodal. Recent parametrizations of the Skyrme force, 
e.g. SLy230a, NRAPR or LNS, showed a similar behav- 
ior We confirm this behavior both within the CP 
and the TF calculation. 

One important conclusion obtained at T = MeV is 
the fact that the transition density for /3-equilibrium mat- 
ter obtained within a TF calculation almost coincides 
with the prediction from the dynamical spinodal. This 
fact should be confirmed at finite temperature. How- 
ever, for symmetric matter or for isospin asymmetries 
not much smaller that yp — 0.3, the TF transition den- 
sity is larger than the prediction of the dynamical spin- 
odal. This proton fraction is of particular interest for 
neutrino trapped matter for which yp ~ 0.3. In this case 
a complete equilibrium calculation should be done. 

The parametrisation GDFM has a very special behav- 
ior with an instability region larger than all the other 
models, for quite symmetric matter. However, for very 
asymmetric matter the instability region is smaller than 



14 



the one of other DDH models and is of the order of the 
NLW models. 

Another important conclusion drawn in the present 
work is the dependence of the pasta phase extension on 
the isospin asymmetry. For the NLWM it is seen clearly 
that the pasta phase extension decreases if the isospin 
asymmetry decreases. For the density dependent mod- 
els a reduction of the proton fraction from 0.5 to 0.3 
almost does not affect the pasta phase or may increase 
it within the CP calculation and the dynamical spinodal 
approach. Within the dynamical spinodal approach this 
behavior is due to the small concavity of the spinodal 
surface for symmetric matter, and the deformation of the 
spinodal due to the presence of protons, electrons and the 
Coulomb field. A smaller fraction of protons contributes 
with less repulsion and gives a larger instability region. 
Of course the presence of electrons shields the proton 
repulsion and therefore the effect is not so strong as it 
would be for charged matter. An adequate description 
of electrons and the Coulomb interaction is important to 
get a correct description of the pasta phase extension. In 
[3| it was shown that the largest pasta phase extension 
occurs when the inclusion of the Coulomb field is done in 
a self-consistent way. 

We have checked which parametrisations fulfill the con- 
straints imposed in fb] for the derivative of the sym- 
metry energy and the transition density and pressure. 
While the density dependent hadronic models are below 
the lower limit for the symmetry energy derivative, the 
NLWM are above the upper limit. However, both NL3 
and NL35 together with GDFM fall within the transi- 
tion pressure/density limits while TW and DDH5 have 
too large transition densities. The GDFM parametriza- 
tion is the one that satisfies the constraints of more 
closely. It seems that the relation between both quanti- 
ties, the slope of the symmetry energy and the transition 
density, is not, in fact, model free. 

If we had considered the limits on the transition den- 
sity defined in [14], 0.25 < Pt < 0.65 MeV/fm^, all mod- 



els studied here have too small transition pressures except 
NL3, that lies just at the lower border, and TW. In Q 
it was shown that a larger slope L gives rise to a smaller 
transition density and transition pressure. This feature is 
seen when models within the same framework are consid- 
ered, namely NL3 and NL3^. However, no clear trend is 
seen among the DDH models. In [23] it was shown that 
the slope of the symmetry energy of the NLW models 
differs from the one of density dependent hadronic mod- 
els. The parametrizations of Skyrme forces, like SLy230a, 
NRAPR or LNS, which have values of L at saturation 
close to the ones of DDH models, also have the transition 
densities close to the ones of DDH models and above the 
limit 0.065 fm^'^ imposed in This seems to show that 
a more complete relation between the transition pres- 
sure and transition density and the equilibrium isovector 
properties of asymmetric nuclear matter have to be ob- 
tained. This could include constraints on the slope and 
compressibility of the symmetry energy at subsaturation 
densities. For instance, the NLW models have positive 
compressibilities of the symmetry energy at subsatura- 
tion densities above 0.05 fm""^ while the DDH models and 
the recent Skyrme parametrizations have negative com- 
pressibilities. Another point which should be referred is 
that the data obtained from isospin diffusion in heavy- 
ion reactions correspond to isospin asymmetries that are 
far from the ones occurring at /3-equilibrium matter. For 
these large asymmetries we expect that the contribution 
from terms beyond the parabolic approximation for the 
isospin dependence of the energy density of nuclear mat- 
ter becomes important [23j . 
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